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Abstract 

We analyze quantum Monte Carlo data in the vicinity of the quantum tran- 
sition between a Neel state and a quantum paramagnet in a two-layer, square 
lattice spin 1/2 Heisenberg antiferromagnet. The real-space correlation func- 
tion and the universal amplitude ratio of the structure factor and the dynamic 
susceptibility show clear evidence of quantum critical behavior at low temper- 
atures. The numerical results are in good quantitative agreement with 1/iV 
calculations for the 0(N) non-linear sigma model. A discrepancy, reported 
earlier, between the critical properties of the antiferromagnet and the sigma 
model is resolved. We also discuss the values of prefactors of the dynamic 
susceptibility and the structure factor in a single layer antiferromagnet at low 
T. 
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The two-layer, square lattice spin 1/2 Heisenberg antiferromagnet has recently attracted 
some attention for a number of reasons. It has been proposed that its properties should 
help explain neutron scattering, photoemission and relaxation rate data in YBCO and 
Bi — 2212 cuprate superconductors^Hl. On the theoretical side, it is among the simplest, 
two-dimensional Heisenberg spin systems which display a quantum phase transition between 
a Neel state and a quantum paramagneti. Quantum Monte Carlo (QMC) simulations on this 
model do not suffer from a sign problem, making it possible to obtain precision data on the 
critical properties of the quantum transition. It therefore constitutes an attractive setting for 
testing our understanding of transitions in two-dimensional Heisenberg spin systems. This is 
particularly important as there are conflicting reports on the observation of quantum-critical 
behavior in a single-layer antiferromagnetS^i. 

A popular approach for the analysis of the collective properties of two-dimensional, square 
lattice antiferromagnets has been to model the long-wavelength, long-time spin fluctuations 
in terms of those of the 0(3) non-linear sigma modetf . The main difference between the 
antiferromagnet and the sigma model is that the microscopic action for the former contains 
Berry phases, associated with the commutation relations between different components of 
the Heisenberg spins. The Berry phases however cancel between the two sublattices for 
smooth spin configurations^. For the two-layer antiferromagnet there should be a further 
cancellation between the two layers, which should occur even for singular saddle points (pairs 
of "hedgehogs" on the two layers) of the effective action. There are, therefore, strong theo- 
retical reasons for believing that the quantum transitions in the two-layer antiferromagnet 
and the sigma model should be in the same universality class. 

In a recent paperS reporting QMC results on the two-layer antiferromagnet near the 
critical point, it was found that while the measurements were qualitatively consistent with 
the predictions of the sigma model, there was a discrepancy in a coefficient in the temper- 
ature (T) dependence of the correlation length (£). The origin of this discrepancy was not 
understood, and it was not clear whether the theoretical picture based on the sigma-model 
should be questioned. On the other hand, high temperature series studies! found good 



agreement with the predictions of the sigma-model at somewhat higher temperatures. 

In this paper, we re-examine this issue by improving the accuracy of the numerical 
calculations, and performing a more detailed analysis of the data. As shown below, we 
now find complete consistency between the measured universal critical properties of the 
two-layer antiferromagnet and the sigma model, and understand the origin of the earlier 
discrepancy^. For comparison, we will also present some additional data on the correlation 
function in a single layer antiferromagnet at low temperatures, when the system is in the 
renormalized-classical regime. 

We start by presenting some new analytical results. Consider the spin correlation func- 
tion in real space, G(r), where (— l) r G(r) is the Fourier transform of the static structure 
factor, 

S(k)=TY,x(k,uj m ) } (1) 

<^777 

where x is the dynamic susceptibility, and u m is a Matsubara frequency. Here and below 
we set ks = ft = 1. Results for x(k,u m ) in the quantum-critical region have been presented 
in Ref 5 — using these results, we obtain at the critical coupling, but T finite 

Here a; is a rescaled momentum, J is a Bessel function, c is the spin-wave velocity, and p s 
and N Q are spin stiffness and sublattice magnetization respectively, measured just away from 
the critical point on the ordered side — the ratio Nq/pI +v remains finite upon approaching 
criticality. Finally, is the universal crossover function, which we have computed in a 
1/N expansion for the O(N) sigma-model for < x < 5 (previously!, this function was 
computed only in the limits of small and large x). We found that to the accuracy needed 
here, S(x) to order 0(1/N) can be modeled by the same functional form as in the N = oo 
result, 



where n x = (exp \/x 2 + m 2 — is a Bose function, but with m « 1.04 instead of the 
N = oo, m = 21og[(v / 5 + l)/2] = 0.9624. The difference between the 1/iV result for H, 
extended to the physical case of iV = 3, and eq. (0) is less than 1% for x < 2 and does 
not exceed 2% for larger x. This form of H indeed implies that the argument of the Bessel 
function can be rewritten as 1.04xr/£ (where £ _1 = mT/c). 

The two terms in (§) yield the same exponential decay of the spin correlations at large 
distances, but with different prefactors. The term without a Bose function is a purely 
quantum contribution, which yields e~ r ^/r, while the integration over x in the second 
term in (§) yields a classical result e r// '»/(r/^) 1 / 2 . In the earlier report, the QMC data near 
the critical point were fitted using only the first quantum contribution to G(r) at large 
distances^. Meanwhile, we have checked that the dominant contribution to S(k) in (|]) 
actually comes from the few first terms in the frequency sum. Consequently, the piece in 
G(r) associated with the classical fluctuations is at least as important as the piece related to 
quantum fluctuations. Moreover, as the QMC data were obtained for finite-size systems, we 
found it useful to extract the correlation length from the fit to the functional form of G(r) 
at moderate distances of 5 — 10 lattice spacings. At such distances, G(r) is not necessarily 
well approximated by the exponential form, and one should use the full eq. (@) for numerical 
comparisons. 

The advantage of introducing the correlation length into the argument of the Bessel 
function is that eq. (0) no longer contains the spin-wave velocity which in a two-layer 
system is not known precisely^. Once £ is obtained from the fit to the QMC data, one 
can use the QMC results for the uniform susceptibility (x u ) an d compute the dimensionless 
ratio QiT) = £ _1 (T)/(Tx u (T)) 1 / 2 . It is predicted^ that at the critical point, Q(T) takes 
a T-independent universal value — using the result for \u we obtain Q w 1.92m = 2.0 (at 
N = oo, Q = 1.64). Away from the critical point, we expect that as T — > 0, Q vanishes 
as ~ exp(— 2ti p s /T) on the Neel ordered side, and Q diverges as ~ exp(A/2T) on the 
quantum-disordered side (A is the spin-gap). 
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For the numerical calculations we have used a variant of the Handscomb QMC tech- 
niqueEl. All results for two coupled layers shown here are for systems with 2 x L x L spins 
with L = 32. For a single layer, we discuss results for systems of size 64 x 64. 

The QMC data000 for G(r) in the two-layer model at a near critical J2/J1 = 2.55 
are presented in Fig. |T] (J 2 and J\ are inter and intralayer exchange integrals, respectively). 
We fitted the data for G{r) to eq. (@) with £ and the prefactor as the fitting parameters. 
We then used the QMC results for the uniform susceptibility^ and obtained Q(T). The 
results, along with results for systems just above (J2/J1 = 2.6) and below {J%j J\ = 2.5) 
the critical point, are shown in Fig [| The results are in good accord with the theoretical 
predictions discussed above. At the critical point J2/J1 = 2.55, Q changes little with 
temperature below T ~ 0.5Ji and remains close to the theoretical prediction, Q = 2.0, (i.e, 
= 1.04T/c) down to the lowest temperature studied. To the best of our knowledge, this 
is the first observation of the universal quantum-critical behavior of the correlation length 
in 2D antiferromagnets. 

Another interesting dimensionless ratio! of the quantum-critical region is W(T) = 
S(n)/Tx(^,0). On the ordered side, as T — > 0, W — > 1 as only the first term in the 
frequency sum in ([I]) is important in this classical regime. At the critical point we however 
have to perform the full frequency summation. Using the 1/N results of Ref 5, we find 
W ~ 1.05m = 1.09. The QMC results for the ratio are shown in Fig. |3|. We see that at 
the transition point, W is T independent and it value is close to the predicted one down 
to the lowest temperatures studied. Notice that earlier studies! have found the same value 
of W in a single-layer antiferromagnet, but at relatively high temperatures, T > O.6J1, 
when nonuniversal corrections to x an d S are relevant, and one can only hope that they are 
canceled out in the ratio S/Tx- 

The prefactor A = Nq (Ji/2T[p s ) (T/2np s ) v extracted from the fits to the correlation 
function is shown in Fig. [|. We see that A w 0.063 is almost T-independent, as it should be 
given that n is very small. This prefactor can also be extracted from the product V(T) = 
S(tt)xu{T) which again does not contain the spin-wave velocity. In the 1/N expansion!, 
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we obtain V(T) = 1.7A. Averaging the QMC data for V at T = 0.3Ji — 0.5Ji, we obtain 
V ~ 0.107 which yields A = 0.063 in complete agreement with the above result. 

Finally, we extracted m and A directly from the fit of the QMC data to eq. (0). We 
used c = I.8J1 which is the average of c = 1.7 Ji and c = 1.9Ji obtained from the fits to 
the data for the uniform susceptibility^ and the ratioS S(k)/Tx{k, 0) at k < ir. Using m 
and A as fitting parameters, we obtained A = 0.063 and m = 1.02; the latter is again very 
close to the theoretical prediction. The A factor was also extracted from QMC calculations 
of the static structure factor S(tt). At the lowest available T ~ 0.3Ji, we obtained using 
c = 1.8 Ji, A = 0.063, in complete agreement with the above results. 

For comparison, we also discuss some new low-temperature data for a single-layer an- 
tiferromagnet, particularly for the prefactor of the correlation function. At low enough 
temperatures, the physics of a single-layer antiferromagnet is governed by thermal, classical 
fluctuations. The N = oo result for H is now qualitatively wrong, as the 1/N corrections 
are logarithmically divergent, and the series for the 1/N terms gives rise to an extra power 
of T in the structure factor. The real space correlator G(r) now has the form 



G(r) = Z 



N 



NqT f {N-2)T \ 1/{N ~ 2) 

ps V 27r Ps / 

'xrT\ 



f°° fxr±\ 

I dx xf(x/m) J f J i (4) 

where f(y) is a universal scaling function (/(0) = 1), and is the overall renormalization 
factor. The behavior of f(y) at large y is rather complex, but we actually need to know 
f(y) only for y = 0(1) as the momentum integration in @ is confined to x ~ m. For such 
y, virtually all 1/N corrections to the mean-field expression for /; f(x/m) = (x 2 + m 2 ) -1 , 
can be absorbed into the renormalization of the mass!, so that for practical purposes we 
again can use the mean- field functional dependence of /. This is consistent with the earlier 
analysis of Tyc et aZ.0, who introduced a semi-phenomenological one-parameter scaling form 
for f(y). For y = 0(1), this scaling form is well approximated by the mean-field expression. 
The prefactor Z^ was obtained^ analytically to first order in 1/N. To estimate Z 3 , we 
assume that the 1/N corrections collapse into a single exponent (this is known to be true 
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1 + 0.188/iV, while for iV = 3, we have Z z ss 2.15, 



3/(JV-2) 



r(i + i/(Ar-2))' 



2 (2-JV) 



(5) 



which is larger than the first-order result but still substantially smaller than the previous 



The QMC data for the real-space correlation function and their fit to the eq. (§) are 
shown in Fig. [|. There is no need to exclude the spin-wave velocity, as in a single-layer anti- 
ferromagnet c and other input parameters are known to high accuracy^!. The mass extracted 
from the fit at low temperatures agrees well with the results of other numerical studies^ and 
with the theoretical formula for the classical regioni'0, though perfect agreement with the 
theory is reached at slightly smaller p s . The prefactor is now A = (A^/27rp g ) [Z 3 (T/27rp s )], 
i.e., it should decrease linearly with T. The A inferred from the fit indeed decreases with 
T, though there are not enough points at low T to obtain the T dependence. The value 
of A at the lowest temperature in the simulations, T = 0.3Ji is A = 0.07, while from 
the theoretical formula we obtain A « 0.05. This indicates that the actual value of Z 3 
is probably about 40% higher than in Another way to find A is to consider the ra- 
tio X = S(n)Jf/(T£) 2 = N$Z 3 Jf/(2irp 2 s ). The QMC result for this ratio is X « 1.34 at 
T = 0.3Ji, while the theory gives X = 0.98, which is again less by about 40%. 

To summarize, in this paper we presented the first numerical evidence for the universal 
quantum-critical behavior of the correlation length in a 2D antiferromagnet near the criti- 
cal point at low temperatures. We compared the results for the universal parameter ratios 
£ -1 1 '(Txu) 1 an d S(ti)/[Tx(tt,0)] with QMC data, and in both cases obtained good agree- 
ment with the theory based on the 0(3) sigma-model. The prefactor A extracted from the 
fit is nearly T independent which agrees with the theoretical prediction A c<. T v . For com- 
parison, we also presented data for the correlation function of a single-layer antiferromagnet 
at low T. We found good agreement with the theoretical result for the correlation length 
in the renormalized-classical region, but the prefactor in the correlation function disagrees 



estimate of Tyc et al. $—Z, w 6.63. 
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with the classical formula by about 40%. 
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FIGURES 

FIG. 1. QMC results 

for the staggered spin correlation function G(r) = (— l) r ((5f r — S'lrX'^To ~~ ^20)) (^ ne indices 
1,2 refer to the two planes) along r = (r, 0), for an interplane coupling J2 = 2.55Ji at tempera- 
tures T = 0.3 — 0.8 Ji (the top curve is for T = 0.3 J\). The solid curves are fits to G(r) + G(L — r), 
with G(r) given by ([|). 

FIG. 2. QMC results for the universal parameter ratio Q = ip 1 /(Txu) 1 ^ 2 f° r the interplane 
couplings J2/J1 = 2.50 (open circles), 2.55 (solid circles), and 2.60 (open squares). The dashed 
line is the prediction for the quantum-critical regime. 

FIG. 3. QMC results for the universal ratio S(tt)/[Tx(t?)] for a single plane (open circles), 
and two coupled planes with J2/J1 = 2.55 (solid circles). The dashed line is the prediction for the 
quantum-critical regime. 

FIG. 4. QMC results for the prefactor, A, in the correlation function for interplane couplings 
J2/J1 = 2.50 (open circles), 2.55 (solid circles), and 2.60 (open squares). 

FIG. 5. QMC results for the spin correlation function C(r) = (— 1Y{S^Sq) along r = (r, 0) for 
a single- layer antiferromagnet at temperatures T = 0.3 — O.8J1 (the top curve is for T = 0.3Ji). 
The solid curves are fits to G{r) + G(L — r), with G(r) given by (H). 
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